Evanescent states in 2D electron systems with spin-orbit interaction and 
spin-dependent transmission through a barrier 
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We find that the total spectrum of electron states in a bounded 2D electron gas with spin-orbit 
interaction contains two types of evanescent states lying in different energy ranges. The first-type 
states fill in a gap, which opens in the band of propagating spin-splitted states if tangential mo- 
mentum is nonzero. They are described by a pure imaginary wavevector. The states of second 
type lie in the forbidden band. They are described by a complex wavevector. These states give rise 
to unusual features of the electron transmission through a lateral potential barrier with spin-orbit 
interaction, such as an oscillatory dependence of the tunneling coefficient on the barrier width and 
electron energy. But of most interest is the spin polarization of an unpolarized incident electron 
flow. Particularly, the transmitted electron current acquires spin polarization even if the distribu- 
tion function of incident electrons is symmetric with respect to the transverse momentum. The 
polarization efficiency is an oscillatory function of the barrier width. Spin filtering is most effective, 
if the Fermi energy is close to the barrier height. 
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I. INTRODUCTION 



The spin-orbit interaction (SOI) in low-dimensional 
structures attracts a great deal of interest since it opens 
up the possibility to manipulate the electron spin in 
nonmagnetic structures using electrical means.™* 2 - In this 
view, semiconductor heterostructures with 2D electrons 
are very promising since the Rashba SOI is effectively 
controlled 3 -* 4 ^ by varying applied bias or gate voltages. 
In recent years, predominant interest was paid to effects 
appearing when the SOI modifies propagating electron 
modes with energy above the conduction band bottom. 
Suffice it to mention the spin-Hall effect jSiLS or spin ma- 
nipulation in strained semiconductors^. In this paper 
we show that interesting effects of the SOI arise also 
when the electron energy is lower than or near to the 
conduction band bottom and evanescent states are in- 
volved. These states determine electron tunneling. They 
are important in 2D structures with laterally inhomoge- 
neous potential landscape. We find that such structures 
can effectively polarize the transmitted electron current. 

3D tunnel structures, in which spin polarization arises 
due to the SOI, were considered in a number of recent 
works. Zakharova et ali-° studied the interband tunnel- 
ing, Voskoboynikov et alii considered a tunnel structures 
with the Rashba SOI at the interfaces. In these cases the 
electron flow acquires a spin polarization if the struc- 
ture is asymmetric. In symmetric tunnel structures the 
spin polarization arises, if the barrier material is non- 
centrosymmetricalii^ The polarization mechanism, pro- 
posed by PereP et a l 12 ' 13 , consists in a spin-dependent 
renormalization of the electron effective mass owing to 
the Dresselhaus SOI in the barrier. However, all these 
structures have a common property restricting their ca- 
pability to generate spin polarization. The polarization is 
absent if the electron current is perpendicular to the bar- 
rier. In other words, for the spin polarization to appear 
the momentum distribution function of incident electrons 



must be asymmetric with respect to the momentum com- 
ponent parallel to the barrier. 

The effective mass renormalization occurs if the Hamil- 
tonian of the SOI is quadratic in longitudinal momen- 
tum. However, the dispersion relation of electrons in 
the presence of the SOI is generally much more compli- 
cated and therefore a more careful analysis of the com- 
plex band structure and evanescent states should be car- 
ried out to study the spin-dependent tunneling. In the 
3D case, such calculations were recently carried out for 
some specific materials and qualitatively new features 
were found. 14 ' 15 ! 16 

2D tunnel structures are scantily studied to date. In 
particular, as far as we know, even the complex band 
structure of 2D electrons was not explored. Though 
the presence of evanescent states is obvious, only a few 
of works touched upon these modes. Usaj, Reynoso 
and Balseiro 17 ' 18 attracted evanescent states to study 
the electron scattering at the edges of 2D samples, but 
the total spectrum of evanescent states was not consid- 
ered. The importance of evanescent modes in quasi- 
one-dimensional systems in the presence of the SOI was 
pointed out in a number of works . 19 ' 20 ' 2 ' 22 

Khodas, Shekter, and FinkePstei n 23 i 24 studied the elec- 
tron beam propagation in 2D electron gas with spatially 
inhomogeneous SOI. They considered the transmission 
through a strip, in which the SOI strength differs from 
that in the rest of the 2D electron gas, to show that an ini- 
tially unpolarized beam splits into two beams with differ- 
ent spin polarizations propagating in different directions. 
The consideration was restricted by propagating states 
since only the case of uniform potential landscape was 
studied. Spin-dependent reflection of electrons from a 
lateral barrier in 2D system was observed in InSb/InAlSb 
heterostuctures 2 ^ and described theoretically in Ref. [26|. 
Silvestrov and Mishchenko 2 ^ demonstrated the possibil- 
ity to achieve spin-polarized currents in a 2D system with 
smooth potential barrier and spatially-uniform SOI by 
considering propagating modes within semiclassical ap- 
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proach. 

In this paper we show that strong polarization effect 
appears in 2D structures when electrons pass through 
a lateral potential barrier, in which the SOI is stronger 
than in the outside 2D electron gas regions (reservoirs). 
The polarization arises even if the electric current is nor- 
mal to the barrier. The highest polarization is attained 
when the electron energy is close to the conduction band 
bottom in the barrier. In this case the fact becomes im- 
portant that the SOI effectively splits the barrier height 
so that some part of electrons passes through the bar- 
rier via propagating states while others do this via the 
evanescent modes. Since the spin and orbital degrees of 
freedom are coupled, rather strong spin filtering occurs. 

We study the total spectrum of electron states in 2D 
bounded system with SOI to find that there are two type 
of evanescent states. The first-type states are character- 
ized by an imaginary longitudinal wavevector. They fill 
in a gap in the propagating state spectrum. The states 
of the second type lie in the forbidden gap. They are de- 
scribed by a complex wavevector. The electron tunneling 
through a lateral barrier with SOI via these evanescent 
states exhibits unusual features, such as an oscillatory 
behavior of the transmission coefficient with the barrier 
width. But of most interest is the spin polarization of 
the electron current. The polarization efficiency is high 
enough even if the distribution function of incident elec- 
trons is symmetric with respect to the transverse mo- 
mentum. We explore the polarization efficiency in a wide 
range of electron energy to find that most effective spin 
filtering occurs if the Fermi energy is close to the barrier 
height. 

The paper is organized as follows. In Sec. |TT] we de- 
scribe the complex band structure and total spectrum 
of electron states. Sec. IIIII is devoted to the tunneling 
through a barrier with SOI. In Sec. IIVI the electric and 
spin currents through the barrier with SOI are considered 
for a wide range of the Fermi energy in the reservoirs. 
The results obtained for the Rashba SOI are generalized 
to the Dresselhaus SOI in Sec. [Vj We summarize and 
conclude in Sec. IVII 



II. COMPLEX BAND STRUCTURE OF A 2D 
ELECTRON SYSTEM WITH SOI 

We start by considering the total spectrum of propa- 
gating and evanescent electron states in 2D electron gas 
with the Rashba SOI. The Hamiltonian is2£ 



where s = +, — stands for spin states, k — (k x ^ ^y)? is 
a constant. The eigenenergy £k,s and the spin function 
X s (k) are defined by following equations: 



(Cm - fc 2 )x«(k) - 2a{k y + ik x ) = 
-2a(k y - ik x ) X s(k) + (Cm - fc 2 ) =0 



(3) 



where o = ma/h 2 is a characteristic wavevector of SOI, 
k 2 = k 2 + k 2 and Cm is the normalized eigenenergy: 



6. 



2m£k,. 



Using Eqs ([3]) one obtains the dispersion equation 

(Cm - k 2 ) 2 - 4a 2 k 2 = (4) 



and the spin function 
Xs(k) 



2a ky 



Cm ~ k 2 



(5) 



Let us analyze the dispersion equation. To be specific, 
assume that the system under consideration is infinite in 
y direction and has a boundary in x direction, i.e. the 
system is semi-infinite or finite in the x direction. In 
this case, the y component of the wavevector k y is real, 
while k x is generally complex, k x = k' x + ik x . Divid- 
ing the real and imaginary parts of Eq. ([4]) we obtain 
an equation set determining the energy ( as a function 
of k' x , k xl k y and trajectories in the (k' x , k' x ) plane, along 
which ((k x ,k x ,k y ) is real: 

(C - k 12 + k" 2 ) 2 = 4a 2 (k' 2 - k" 2 ) + Ak' 2 k" 2 , (6) 
(C - k 12 + k" 2 )k'k" = -2a 2 k'k" , (7) 



k 2 y . 



where k' + ik" = ^J(k' x + ik'J. 

Eq. ([7]) possesses a solution in the following cases: 
i) k' = 0, ii) k" = and iii) ( - k' 2 + k" 2 = -2a 2 
for k', k" ^ 0. The first case contradicts to Eq. ([6]) and 
hence cannot be realized. The second case generates two 
branches of the solution, which can be found after divid- 
ing the real and imaginary parts of Eq. ^ : 

1. k" = 0, 



> 2 

Cm = - a 



a ± \/k 2 



X.(k) = ± ky + lK 



(8) 



1,2 i hi 2 

K y + K x 



HR = —(pl+p 2 y ) + ^{p y a x -p x a y ), (1) 2. k' x = 0, |A£'| < \k y \, 



where a is the SOI constant, a x ,a y are Pauli matrices. 
The eigcnfunctions are 



(2) 



a±\lkl - k'l' 



k — k" 

x.(k) = ± ky K 



1,2 _ I,// 2 

K x 



(9) 
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spin function are 




FIG. 1: (color online). Total spectrum of 2D electron gas 
with SOI. Solid lines 1+ and 1— are spin-splitted propagat- 
ing modes of branch 1; lines 2+ and 2— are evanescent states, 
branch 2; lines 3L, 3R are evanescent states in the forbidden 
gap, branch 3. Dashed lines 1 , 2' and 31/ , 37?' are the trajec- 
tories corresponding to these branches on the complex plane 
\h x , h x ) . 



The third case gives one further branch: 

3. This branch is defined for k' x ,k x belonging to the 
folowing trajectory 

k' x 2 kf + a\kl + K 2 - kf) a 4 = (10) 
in the complex plane (k' x ,k x ). The eigenenergy and the 



X.0O = 



y 2 h" 2 



ky k^, ~\~ ^kj, 
1 a 2 + ik' x k'J. 



(11) 



The total spectrum is shown schematically in Fig. Q] 
where the energy £k, s is presented as a function of k' x , k" 
for given transverse momentum k y . The form of all three 
branches is different in the cases \k y \ < a and \k y \ > a. 

Branch 1 describes the propagating states with real k x . 
This branch is splitted by spin. The energy gap, which 
opens at k x — for k y ^ Q, depends on \k y \. 

Branch 2 describes purely decaying evanescent states 
defined on the imaginary k x axis. This branch con- 
nects the propagating state branches along the imaginary 
axis. The spectrum of branch 2 is also splitted by spin. 
The energy minimum s m — —ha 2 /(2m) is common for 
branches 1 and 2. The minimum is attained on branch 
1, if \k y \ < a, or on branch 2 in the opposed case. 

Branch 3 corresponds to evanescent states in the for- 
bidden gap, £k,s < £m- They are described by a com- 
plex longitudinal wavevector and therefore can be named 
"oscillating" evanescent states. The trajectories, along 
which these states are defined, obey the following equa- 
tion: 



(a 2 -k' x 2 )(a 2 + kf) 



i 2 e, . 



(12) 



They are shown in Fig. [5] There are trajectories of two 
types. If \k y \ < a, the trajectories (lines 1, 2) intersect 
the imaginary axis k x . In the vicinity of the intersection 
point, \k x \ ^ \k' x \ and hence the wavefunction oscillates 
with the distance faster than decreases. If \k y \ > a the 
trajectories (lines 4-6) intersect the k' x axis. 



k'x/ak 




FIG. 2: (color online). Real-energy trajectories along 
which the "oscillating" evanescent states are defined. k y 
is fixed for each line: lines 1-6 correspond to k y /a = 
0.4,0.8,1.0,2.0,5.0, 10.0. 



It is seen that there are four states for any energy, in 
accordance with the number of the degrees of freedom 
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of the system. At a given energy all the states are dis- 
tinguished by the wavevector components k' x , k" and the 
spin function. The propagating states and the evanes- 
cent states of branch 2 have two subbranches divided by 
the energy, while the evanescent states of branch 3 are 
splitted in the complex momentum plane. 

Though the wavefunctions of the third branch are com- 
plex, they do not carry the current. Using the Hamil- 
tonian |T]) one obtains the following expression for the 
particle current in the state 



4>s 



2m 



(13) 



(14) 



The wavefunctions of branch 3, calculated with using 
Eqs and (jTTJ) , are easily seen to turn the x component 
of the current (|14|) to zero while y component is nonzero. 



III. 



TUNNELING CURRENTS THROUGH A 
BARRIER WITH SOI 



In this section we study the electron tunneling via the 
"oscillating" evanescent states of branch 3. Consider a 
2D electron gas without SOI divided into two semiplanes 
(reservoirs) by a rectangular barrier, in which the SOI is 
present. The barrier height is U and the width is d. 

Let us calculate the transmission probability for elec- 
trons incident on the barrier from the left reservoir. We 
use the presentation of the wavefunctions in the reser- 
voirs in the basis of eigenstates \s) of the cr z matrix. The 
basis functions are \k x ,k y ,s), where k x and k y are the 
wavevector components in the reservoirs, s =1, j. 

In the left reservoir (x < 0) the wavefunction is 



) — \k Xl ky, s) + y j t s s i | k x , ky , s ) , (15) 



where \k x ,k y , s) is the state vector of incident electrons. 
The wavefunction of transmitted electrons (x > d) is 

l^fi, s >=EW#*A,s'>- (16) 



Here r S)S > and t s-s > are reflection and transmission matri- 
ces. 

The wavefunction in the barrier is expanded in the 
eigenstates @ of the Hamiltonian |T]) 



\4 B X,s)= E K tr ,\rK' x ,r'K':,k y 



(17) 



where K x = K' x + iK" denotes the complex wavevec- 
tor in the barrier; r,r' = ± are indexes labeling all four 



evanescent states in the barrier. They are described by 
the spinors in which the real and imaginary parts of 
K x should be taken with different signs. 

The matrices r SjS /, t SjS ' and b s rr , are determined by an 
equation set, which follows from the boundary conditions 
at the interfaces of the barrier and 2D electron reservoirs. 

Boundary conditions for wavefunctions are obtained in 
a standard way by integrating the Schrodinger equation 
over the infinitesimal vicinity of the boundary. These 
conditions are well known for a boundary between re- 
gions with different strength of SO I 23 ' 29 . In our case it is 
necessary to take into account that the lateral potential 
step at the boundary also contributes to the SOI. In the 
transition region, where the potential U(x) varies with 
x, the following additional term should be added to the 
Hamiltonian fl} 



jdU 
h dx 



(18) 



where 7 is SOI constant connected with a by a = —e^Fz, 
F z being the electric field perpendicular to the 2D layer. 
This term having been integrated over the transition re- 
gion gives a finite contribution to the boundary condi- 
tions. Finally one obtains the following equations for the 
spinor (|13p components: 



Vi|+o = <H^=0; 



1 



m(x) 
1 



m(x) |_ dx 
where the parameter 



dip 

-g— ± I3(x)k y ipi ~ a(x)ip 2 
dip 

=F (3(x)kyip 2 + a(x)ipi 



+0 

-0 
+0 



= 



0. 



R 2U 



(19) 
(20) 

(21) 
(22) 



describes the Rashba SOI caused by the in-plane field. 
The upper and lower signs in Eqs (f2"0"|) , (|2"Tj) correspond 
to the boundaries at which U(x) increases or decreases 
with x. 

Applying these boundary conditions to the system un- 
der consideration we put a = in the 2D electron reser- 
voirs and keep a ^ in the barrier. For simplicity, we 
ignore the difference in the effective masses of electrons 
in the barrier and reservoirs. 

In addition, the wavevectors in the barrier and the 
reservoirs should be matched. The tangential compo- 
nents k v must be equal. The relation of the normal com- 
ponents is determined by equaling the energy E of an 
incident electron to the electron energy in the barrier 



E(k x ,k v ) =U + s(K x ,k y ,s) 



(23) 



where e(K x , k y ,s) is defined by Eqs (|8|) . (|9|) . (fTT]) m accor- 
dance with the energy spectrum branch (1, 2 or 3), which 
is considered. In the present section we restrict ourselves 
by branch 3. 
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Finally one obtains two equation sets for the cases of 
the incident spin directed up and down. Each equation 
set contains eight equations. Dropping elementary cal- 
culations and combersome expressions for matrices t S)S ', 
r StS > and b* r ,, we turn directly to main results. 

First, note that the t SjS < matrix obeys the following 
symmetry relations: 



The spin density of transmitted electrons is 



tfi(K x , k y ) =t u (K x , -k y ) , 
tn(Kx, k y ) = - t u (K x , -k y ) . 

The particle current through the barrier is 

jk x ,k y ,a — T s (k x , k y ) , 



(24) 



(25) 



where T s is the transmission coefficient, T s = Y] , |i S! s'| 2 , 
which has a meaning of the probability for an electron to 
tunnel with any spin in the final state. The peculiarity of 
the tunneling in the presence of the SOI consists in the 
involvement of four interfering states with different spin 
structure. 

One of unusual consequences of this fact is an oscilla- 
tory dependence of the tunneling coefficient on the bar- 
rier width. This feature is demonstrated in Fig. [3] The 
oscillations exist when the electron energy is close to the 
top of the third branch: U — E so — E <C E so , where 
E so = h 2 a 2 /(2m) is a characteristic energy scale of SOL 
The oscillations fade away as the energy decreases deep 
into the forbidden band, because K" exceeds K' x . The 
oscillations disappear also if the tangential momentum 
k y goes to zero. 




Pk x ,k y 



i5> 



(R) 



(26) 



With using Eq. (fT6|) . the spin density takes the form: 



3yS )t s s /tg — g' -j- G Z S \t S 



where |C| 2 is a normalization constant. 

The components of the spin polarization are shown in 
Fig. H] as functions of the transverse momentum k y for a 
given energy of incident electrons. The main property of 
the acquired polarization is that the x and z components 
of the spin polarization are odd functions of k y , while the 
y component is an even function of k y . This property is 
independent of the electron energy, the barrier height, the 
SOI strength, and thus is universal for the Rashba SOI. 
Two consequences follow from this fact. First, the un- 
polarized electron follow acquires spin polarization even 
if the distribution function of the incident electrons is 
symmetric with respect to the tangential momentum di- 
rection, the polarization being directed parallel to the 
barrier. Second, if the distribution function is not sym- 
metric about k y , the spin polarization arises also in the 
x and z directions. 
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FIG. 4: (color online). Spin polarization of transmitted elec- 
trons, £l = 2fj, XtV ^/h\C\ 2 . The curves are marked by letters 
corresponding to the polarization components. The incident 
electron energy E = 9.5 E so , the barrier height U = 11 E so , 
the barrier width d — a, the lateral SOI parameter j3 = 0.1. 



FIG. 3: (color online). Tunneling coefficients for incident 
electrons with spin up and down as functions of the barrier 
width. U = 2E ao , E = 0.95£ SO , k y = 0.49a and (3 = 0.1. 

Fig. [3] shows clearly that the barrier filters incident 
electrons by spin. This process depends evidently on 
the incident angle and the energy of electrons. Consider 
the spin polarization of transmitted electrons in more de- 
tail. Let an unpolarized electron flow with the wavevector 
(k x ,k y ) is incident on the barrier. The flow consists of 
two components \k x ,k y ,1) and 1^,^,1) with opposed 
spins. The transmitted flow acquires spin polarization. 



Detail analysis shows that the spin polarization is 
caused mainly by the SOI in the barrier, as it is described 
by the Hamiltonian ([T]). The SOI at the boundaries of 
the barrier (described by Eq. (|18p) does not essentially 
affect the results if the parameter (3 defined by Eq. (|22|) 
is small. This case corresponds to realistic situation in 
experiments. Numerical estimations for InAs quantum 
well (a - 6 x 1(T 9 eVcm, F z ~ 10 5 V/cm, U ~ 20 meV) 
give j3 ~ 0.1. If j3 > 1, the boundary SOI changes the 
value of polarization, but main features (such as the sym- 
metry relations, the oscillatory behavior of the tunneling 
coefficients) remain qualitatively similar. 
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IV. SPIN POLARIZATION OF ELECTRONS BY 
A BARRIER WITH SOI 

In this section we turn from the separate electron states 
to the total spin polarization produced by all electron 
states contributing to the current through a barrier with 
SOI. In addition, we consider a wide range of the inci- 
dent electron energy to include all three branches of the 
electron spectrum. This allows one to find out condi- 
tions under which the electron current is polarized most 
effectively. 

Consider an electron current directed normally to the 
barrier. For definitcncss, let the current be caused by a 
voltage V applied across the 2D electron reservoirs, as it 
is shown in Fig. [5l The incident electron states, which 
contribute to the current, occupy an energy layer near 
the Fermi level. They are located in a semi-ring in the 
k x ,k y plane, shown in the insertion. The particle and 
spin currents through the barrier are determined by sum- 
ming over all these states. We carry out this calculation 
for various positions of the Fermi level Ep relative to the 
barrier height U to find the spin-polarization efficiency 
as a function of Ep. 




FIG. 5: (color online). A schematic view of the barrier with 
SOI. Lines L and R show the electron dispersion relations in 
the left and right reservoirs. Lines 1,2,3 image schematically 
corresponding branches of the dispersion relation in the bar- 
rier. The insertion is the k x ,k y space, in which the shaded 
region shows the electron states contributing to the current. 



For simplicity suppose that the voltage is small com- 
pared to other energies, eV <C U, Ep, E so . This simplifi- 
cation allows one to restrict the summation by the inte- 
gration in k space over the azimuthal angle at a given en- 
ergy E. The integration is conveniently to carry out over 
k v , but in doing this the fact should be taken into account 
that the set of branches, which must be used at given k y 
and E, can change with varying k y . This happens be- 
cause the gap Aei2 between spin-splitted subbranches of 
propagating modes (curves 1+ and 1- in Fig.[T|) increases 
with k y (Ae 12 = 2h 2 a\k y \/m) and furthermore the form 
of the dispersion curves describing the propagating (1+, 



1-) and evanescent (2+, 2-) modes changes. The switches 
between the actual branches occur when any extremum of 
the dispersion curves (which is a function of k v ) coinsides 
with the energy E. Physically this means that electrons 
incident on the barrier at different angles and in different 
spin states feel different effective barrier height. 

Fig. [S] presents a diagram showing which branches of 
the dispersion relation are accessible for electrons with 
given E and k y . Within each region bounded by thick 
lines, there are two different branches with two funda- 
mental eigenfunctions on each ones or one branch with 
four solutions. 





\k y \/a 

5 






2 ~ ,2 ~ A 




3R,3L 


s 3 


l-,2+ 




' 1-1- , 


C 






— " 



-1 \ 1 772) o, w 1 
^lAJ (E F -U)/E S0 



FIG. 6: The diagram of the spectrum branches accessible 
for electrons with energy E and tangential momentum k y . 
Thick lines divide the (E,k y ) plane into 6 regions, in which 
corresponding branches are denoted by numbers 1,2,3 and 
chirality indexes + , — in accordance with Fig. \T\ Line L is 
defined by equation k y /a = {^(E - U - E ao )/E ao + 1 ± l) 2 



For each region of the diagram the wavefunctions are 
determined in the same way as in Sec. IIIII The only dif- 
ference is that the eigenfunctions of the branches spec- 
ified in the diagram are to be used in Eq. (|17p instead 
of the wavefunctions of branch 3. As a result of these 
calculations the transmission t s ^ s >(E,k y ) and reflection 
r StS '(E,ky) matrices are obtained for the whole (E,k y ) 
plane. Using t StS i(E,k y ) we find the particle and spin 
currents in the right electron reservoir. 

The particle current is found by the summation of 
Eq. (|2"5j) over all states of incident electrons: 



eV 
2irh 



dMl*TT| 2 + l*u| 2 + |ilT| 2 + l*u| 2 ) . (28) 



where the integration symbol implies also the summation 
over all regions of the (E, k y ) plane which fall within the 
(—kp,kp) interval, with kp being the Fermi wavevector 
in the 2D reservoirs. 

The transmitted spin current in the left reservoir is 
defined in a standard wa y 30 ' 31 using the following ex- 
pression for the current in a state \k x ,k y ): 



1 



(29) 
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where i = {x,y} denotes the velocity components in the 
plane, j = {x, y, z) denotes the spin components in 3D 
space. In the case under consideration, the x component 
of the total spin current is 




where three lines in the RHS correspond to the x, y, z 
components of the spin polarization for the spin current 
directed along x axis. 

The efficiency of spin polarization is characterized by 
the ratio of the spin current to the particle current 



T>j = -J • (31) 

Using the symmetry relations ([24]) and Eq. ([30]) one 
finds that x and z components of the spin current are 
absent in the case of the Rashba SOI, and only the y 
component is nonzero. Of course, this is a consequence 
of the symmetry of the distribution function with respect 
to the sign of k y . If the current had not be perpendicu- 
lar to the barrier, the polarization would appear also in 
the x and z directions. The polarization efficiency turns 
out to be sufficiently high. The dependencies of the po- 
larization on the Fermi energy and the barrier width are 
nontrivial because they reflect a complex structure of the 
electron spectrum. Below two most significant results are 
considered. 

The dependence of V y on the barrier width is shown 
in Fig. [7] for energies below the top of the energy band 
where the evanescent states of branch 3 exist in the bar- 
rier, Ep < U — E so = Eh- The polarization efficiency is 
seen to oscillate with d because of the interference of four 
oscillating evanescent modes. The effect is rather strong 
and becomes the stronger the closer is the energy to the 
band top. The oscillation period is of the order of n/a. 



V 




FIG. 7: Dependence of the polarization efficiency on the bar- 
rier width for the Fermi energy below Eh for tunneling through 
branch 3. U = 5E so , E F = 3.5E SOy f3 = 0.1. 

The dependence of V y on the Fermi energy is presented 
in Fig. [5] for a wide energy range including the energy 
both below and above the barrier. It is seen that there 



is a critical energy E = 8E so , which coincides with the 
top energy of the third branch band, Ef, = U — E so , in 
the barrier. At energies well below E U1 the polarization 
efficiency V v increases with the energy. In the vicinity 
of the threshold E u an oscillatory behavior appears as a 
result of the interference of four slowly decaying waves. 
Above this energy the transmission process goes via the 
propagating states of branch 1 and the evanescent states 
of branch 2. At Ep slightly higher than E u the polar- 
ization efficiency attains the highest value. With further 
increasing the energy the efficiency V y decreases. It is 
worth noting that the highest efficiency of the spin po- 
larization slightly depends on the SOI strength, but the 
barrier width, at which this high polarization is attained, 
varies inversely with the SOI constant. 

Fig. [5] demonstrates also how V y changes with the bar- 
rier width. Increasing d above ~ n/a leads to more 
pronounced oscillatory dependence of V y on the energy 
caused by the four-wave interference. 



0.6 n 




Fermi energy, EJE 

FIG. 8: (color online). Dependence of the polarization effi- 
ciency on the Fermi energy in 2D electron reservoirs. U = 
9E so for d = 5/ a (solid line 1) and d — 3/ a (dashed line 
2); (3 — 0.1. Vertical dashed line marks the threshold energy 
Eh = U — E so , which divides the evanescent states of branch 
3 (on the left of E B ) from the states of branches 2 and 1. 

The physical mechanism, owing to which the polar- 
ization of normally incident electron current appears, is 
connected with the splitting of electron waves in the bar- 
rier because of the SOI. Let us consider first a simplified 
case of semiinfinite barrier region with SOI. Let the en- 
ergy is high enough so that electrons occupy propagating 
states (branches 1+,1-). The incident electron states can 
be represented in terms of two chiral modes: 

= / ±x\ = k y + ik x 

V 1 / ' ^jkl + ki 

In the barrier region each electron beam splits into two 
beams, which propagate at different angles and have dif- 
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ferent chiralities: 



(i) 

A. 1,2 



Xl,2 
1 



Xl,2 



± 



■iKi 



h 2 



K" 2 

-"■1.2 



(32) 



where K\ is ^-component of the wavevector of the upper 
mode with positive chirality xi and K 2 corresponds to 
the lower mode with negative chirality \2- It is essential 
that K\ < K 2 . The incident and refracted beams as 
well as their spin polarizations are illustrated in Fig. [9] 
The transmitted beam amplitudes are A+ and B+ for the 
incident beam with positive chirality, and A_ and B_ for 
negative chirality. The x and y components of the spin 
polarization in the barrier are estimated as 

S x cx (|A + | 2 + |A_| 2 )Rc X i-(| J B + | 2 + |B_| 2 )Rc X2 +..., 
S y cx (|A + | 2 + |^_| 2 )Imxi-(|5 + | 2 + |5_| 2 )Im X 2+... , 



At small k y , the spin polarization appears. Since the 
amplitudes are functions of fc 2 , they remain unchanged 
in the first approximation, so that the resulting spin is 
determined by the difference (xi — X2) and 



3 n 



h 2 



h 2 

Ky 



Since K\ < K2, the ^-component of the spin density is 
negative, i.e. the sign of the spin polarization is deter- 
mined by the lower-energy branch of propagating mode. 

If the barrier width is finite, there are four modes in the 
barrier, but the above property remains unchanged. It is 
also kept for other branches of the electron spectrum. 



THE CASE OF DRESSELHAUS SOI 



where the dots denote spatially dependent terms origi- 
nating from the interference. 



B> B - 




2D electrons without SOI 



barrier with SOI 



FIG. 9: (color online). Refraction and spin polarization of 
electron beams incident on the semiinfinite barrier with SOI. 
Solid arrows show the spin polarization in the case of posi- 
tive chirality of the incident beam. Dashed arrows are the 
polarization for negative incident chirality. A± and B± are 
amplitudes for incident beam with k y > 0. The primed letters 
denote the same amplitudes for k y > 0. 

The total spin polarization is determined by the sum 
over all incident angles. To estimate this sum, consider 
the dependence of the spin components on k y . It is clear 
that |^4±| 2 and \B±\ 2 are even functions of k y . Real and 
imaginary parts of \i,2 are seen from Eq. (|32|1 to be cor- 
respondingly odd and even functions of k y . Therefore S x 
is an even function of k y and S y is an odd function, as 
it is illustrated in Fig. [9] Thus the total S x component 
vanishes while 5*^ is nonzero. 

It is easy to find the direction of the total spin for 
the energy close to the barrier height. If k y = 0, the 
amplitudes \A + \ and \B_\ are equal, = \B + \ = 

and Xl = X2 — i- Hence, the total spin density is zero 
because the spins of opposed chiralities cancel each other. 



All the above results are generalized to the case of 
Dresselhaus SOI. For a 2D system oriented along [001] 
crystallographic direction, the SOI Hamiltonian is2£ 



h . 9 2 a 

H D = 7^(Px + Py) + jKVxOx - PyCTy) ■ 



It is well known that Dresselhaus and Rashba Hamilto- 
nians are unitary equivalent^ An unitary matrix 



U 



i 

1 



(33) 



transforms the Rashba Hamiltonian |T]) to the Dressel- 
haus one: Hr = U + HrU = Hd- Therefore, the Dressel- 
haus SOI case does not require separate calculations. It 
is enough to carry out this transformation. Taking into 
account that matrix (|33| transforms the spin matrices as 
follows: a x = —cry, a y = —cr x , and a z — — cr z , we ar- 
rive at the conclusions: (i) the dispersion equation is the 
same as in the Rashba case 

(ii) the spin functions differ from those defined by 
Eq. §5§ by a simple substitution k x ^ k y ; 

(hi) the spin components of transmitted electrons in 
the incident states \k x ,k y ,1) and \k x ,k y ,l) differ from 
those of Sec lIIII by replacements: \i x — > — fly, fi y — > — fx x 
and fx z -> -/i z ; 

(iv) the spin polarization of the current transmitted 
through a barrier with SOI has only x component, if the 
distribution function is even with respect to the trans- 
verse momentum. In particular, Figs [71 [5] are valid for 
the spin polarization normal to the barrier, V x - 

Of course, if the Rashba and Dresselhaus mechanisms 
act simultaneously the results change qualitatively. 



VI. CONCLUSION 

We have found the total spectrum of electron states 
in a bounded 2D electron gas with SOI. It addition to 
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well known propagating states it contains two branches 
of evanescent states. Their wavefunctions decay with the 
distance in the direction x perpendicular to the bound- 
ary. One branch ("purely decaying" evanescent mode) 
is described by an imaginary wavevector. The energy of 
this state is splitted by spin, so that there are two sub- 
branches. They fill in the gap, which opens in the prop- 
agating state spectrum at k y =^ 0. Other branch ("oscil- 
lating" evanescent mode) is characterized by a complex 
wavevector K x . These states lie in the forbidden gap. 

We have studied the electron transmission through 
a lateral potential barrier with the SOI. In the energy 
range, where electrons tunnel via the oscillating evanes- 
cent states, the tunneling reveals unusual features, such 
as an oscillatory dependence of the transmission coeffi- 
cients on the barrier width and the energy. But of most 
importance is the spin polarization of the electron cur- 
rent. The value and direction of the polarization depend 
on the angle of incidence and the energy of incident elec- 
trons. The polarization appears even if the distribution 
function of incident electrons is symmetric with respect 
to the transverse momentum. In this case the polariza- 



tion is directed parallel to the barrier (in the Rashba 
SOI case) or perpendicular to it (for Dresselhaus SOI). 
The highest efficiency of the spin polarization is attained 
when the Fermi energy is close to the barrier height. In 
this case, electrons pass through the barrier partially via 
the propagating states and partially via the purely de- 
caying evanescent states. Under this condition the most 
effective spin filtering occurs. The maximal polarization 
efficiency depends on the barrier height and can exceed 
0.5 if the barrier width is on the order of n/a. 
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